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We find a unitary operator wiiicii asymptotically diagonalizes the Tomonaga-Luttinger hamilto- 
nian of one- dimensional spinless electrons. The operator performs a Bogoliubov rotation in the space 
of electron-hole pairs. If bare interaction of the physical electrons is sufficiently small this transfor- 
mation maps the original Tomonaga-Luttinger system on a system of free fermionic quasiparticles. 
Our representation is useful when the electron dispersion deviates from linear form. For such sit- 
uation we obtain non-perturbative results for the electron gas free energy and the density-density 
propagator. 



I. INTRODUCTION 

Bosonization is a standard approach to the problem of interacting electrons in one dimension (ID) [1, 2]. Bosoniza- 
tion maps the low-energy spectrum of the Tomonaga-Luttinger (TL) model on the spectrum of free bosons. 

In this paper we discuss a new method of treating such system. We explicitly construct a unitary operator U diago- 
nalizing TL hamiltonian. Our unitary transformation maps the TL model on a system of free fermionic quasiparticles. 
The description of the TL model in terms of these quasiparticles has several advantages over bosonization, as we will 
see below. 

The diagonalizing procedure is closely related to the bosonization. The operator U can be thought of as a Bo- 
goliubov transformation in the space of particle-hole pair excitations. Alternatively, one can describe the action of 
[/ as a sequence of bosonization, Bogoliubov transformation which diagonalizes the bosonic hamiltonian and further 
refermionization (fig.l). 

When the dispersion of the physical electrons is linear both the quasiparticles and the bosons offer equally good 
description of the TL model. If the non-linear terms are substantial the free boson representation breaks down. The 
fermionic quasiparticles show more resilience toward deviations from the linear dispersion. They remain free as long 
as the bare interaction constant is sufficiently small. 

This latter property of the quasiparticle representation allows for non-perturbative calculation of the free energy 
and the density-density correlation function for the TL model with the non-linear dispersion. We believe that these 
two results are new. 

When the bare electron dispersion is linear the quasiparticle representation could be used to determine the single- 
electron Green's function. It coincides with the Green's function obtained by other methods. 

The existence of the quasiparticles does not contradict to the fact that the single-electron Green's function of TL 
model has no pole. The quasiparticles in TL model has zero overlap with the physical electrons: \/Z = 0. Therefore, 
it is convenient to think about TL model as Z = Fermi liquid [3]. 

The paper is organized as follows. We diagonalize the TL hamiltonian in Sec. II. In Sec. Ill we offer an intuitive 
explanation to the method. In this section all technical complications are disregarded in order to create an integral 
view. The density-density propagator is derived in Sec. IV. The single-particle Green's function is calculated in Sec. 
V. Sec. VI is reserved for discussion. Certain technical details can be found in Appendices. 

II. DIAGONALIZATION OF TL HAMILTONIAN 

The TL model hamiltonian is given by: 

H = H\^in + Hint, (1) 

i?ki„ = WF / dx(: ^VV'l: -: V'^Vi^j^: , (2) 

J-L/2 ^ ' 

Hint = J dxdx'g{x - x')pi,{x)pn{x'), (3) 
PL,R =: V'l,rV'l,r: ■ (4) 
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The chirality label 'L' corresponds to left-moving electrons, the label 'R' corresponds to right-moving electrons. The 
interaction of the electrons of the same chirality is ignored since up to irrelevant operators such interaction simply 
renormalizes the value of the Fermi velocity v-p. The symbol : . . . : denotes normal ordering of the field operators i/j. 
A brief discussion of the normal ordering procedure we use in this paper is given in Appendix A. 

It is assumed that the cut-off of (1) is infinite. To remove ultraviolet divergences of the theory without cut- 
off we replace usual zero-range interaction (^(a;) = goS{x)) by interaction acting over a finite range. Specifically, 
^(a;) = go^Aix) where Sa{x) is a broadened version of the delta-function: its Fourier transform 5a(9) is such that 
^a{(1) = 1 for \q\ smaller than some quantity A, and SA{q) vanishes quickly for \q\ > A. The parameter A thus defined 
regularizes ultraviolet divergences of our theory. 

It is easy to demonstrate (see, for example, [1]) that the following commutation relations are obeyed: 

[Ppq,Pp'-q'] = 6pp'6gq'pnq (5) 



Op«.V';'(a;)J = Ve-'^^V'M^)' (6) 
where p = +1 for left-moving electrons, p = — 1 for right-moving electrons and 

Ppq = J dxe-"^^pp{x), = (7) 

Using this commutation relations we will show that the unitary operator 

U = e^, (8) 

^ = J2Y1 aq^pp-qp-pq (9) 

diagonalizes the hamiltonian (1) except for the zero mode part. Fortunately, the zero modes 

Np = Ppq (10) 

are decoupled from other degrees of freedom. Also, their contribution to the low-energy spectrum is 0{l/L). 

Since p's are quadratic in tp the above operator U is quartic in fermionic operators ip. In general, it is impossible to 
work with such a non-linear object. In our situation, however, the simplicity of commutation rules (5) and (6) allows 
us to diagonalize the hamiltonian with the help of U. 

In order to transform the interaction term (3) with the operator U it is enough to observe that the action of U on 
the density operator ppq, g ^ 0, is a Bogoliubov rotation: 

UppqU'^ = U{aq)ppq + v{aq)p-pq, (11) 

u{aq) = Uq = cosh(ag), v{aq) = Vq — sinh(aq), (12) 

ul - vl = 1. (13) 

This result is a simple consequence of the commutation relation (5). This identity can be derived by a variety of 
methods. For example, one can differentiate its left-hand side with respect to aq for both p = L and p = R and solve 
the resultant differential equation system. 

The easiest way to transform the kinetic energy term is to notice that the kinetic energy density can be expressed 
as a product of two density operators. The derivation goes as follows. First, we normal order the product of two 
density operators: 

Pp{x)pp{y) = :i)lix)i;p{x)i;l{y)i;p{y):+Sp{x-y):tpp{x)tpl{y):+ (14) 
Sp{x - y): tpl{x)i)p{y): +bp{x - y), 

bp{x) = {sp{x)f , (16) 

As it is explained in Appendix A the normal ordering is used here to isolate explicitly singular terms of the field 
operator products. Now we expand the above identity into Laurent series with respect to powers of {x — y): 

Ppi.x)pp{y) = &p(x - y) -I- —:V'p(a;)V?/'p(a;): -I-—: ■!/'p(a;)VV'p(as-): +(irrelevant operators). (17) 
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In this expansion an irrelevant operator can be recognized by a factor of {x — y)" where n > 0. For example, the first 
term in (14) is an irrelevant operator: 

■.i,l{x)i,^{x)i^l{y)i>p{yy. ^{x-yf:^l{x)i>^{x)Vi>l{x)Vi>^{x):. (18) 

Indeed, by simple power counting one can verify that its scaling dimension is equal to d = 4 < 2. Sending y — > a; in 
eq. (17), we establish the identity: 

^ (: V^(a;)(VVp(a;)): -: {V^l{x))^^{x):) = lim{pp{x)pp{y) - bp{x - y)} . (19) 

Thus, it is permissible to write the hamiltonian (1) in the form: 



H = [PpqPp-q + ^^PpqP-p-q 1 , (20) 

pq ^ ^ ^ 

where is the Fourier transform of g{x). For convenience we explicitly show the zero mode part of the above 
expression: 

Such splitting is useful since our transformation U does not act on the zero modes. We now apply U and choose 
parameters a, in such a way that the term PlPr vanishes: 



9q 

uciiiiiz,u;q — 7 

1 



tanh2a<, = --S^, (22) 
27rWF 



»^-o|l+ / ^ , h ^.' = 2 / ^ J ' ^«5,<0. (23) 

^Jl-{gq/2nVFf ) ^\ yJl-{g,/2nVFf 



The transformed hamiltonian is: 



UHU^ = ^J:J:U+vI + -^u,v,) p,,p,., + ^{Nl + N^) + 'j-N^N^ (24) 

g^O p ^ F / 

= ^ E (-^ + + ;^-.-.) + ^^^^^ (^L + Nl) + |7V,iVH, (25) 

Prom now on we will ignore the zero mode contribution to the hamiltonian: the above formulas clearly show that it 
is small (~ 1/L) for a macroscopic system. Ultimately, we have: 

This expression can be re- written in terms of fields "0. Let us show how this is done in a technically trivial case of 
A ^ oo (which is equivalent to g(x) = go5{x)). When A is infinite the functions Uq, Vq and gq are all constants 
independent of q. Their values are given by (23) with gq = go. Therefore, hamiltonian (27) is equal to 

UHU^ = !^^p^^p^_^, (28) 

pq 

where vp is given by (26). Finally, inverting (19), we obtain 

UHU^ = ivF j dx^p-.tplVipp-. +const. (29) 
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These transformations, apart from minor differences, are equivalent to the refermionization as it is done in [4]. 

In certain situations, however, it is convenient to have finite A. Then a subtler reasoning is required. It is possible 
to prove that: 

UHU'' = ivp / dx ^^Vijjj,: + (irrelevant operators). (30) 

Such proof is rather technical. It can be found in Appendix B. 

Let us agree that operators with the tilde are transformations of corresponding operators without the tilde, that is 
O = WOU. Then the hamiltonian can be expressed in terms of operators -ip as follows: 

H = ivF / dx^^p: V'^V'^^: + (irrelevant operators). (31) 

p 

The above equation is the desired mapping of TL model on the model of quasiparticles ip whose interactions are small 
irrelevant operators. Due to the irrelevance of the interaction the low-energy properties of V' are those of free fermions. 

Using the developed framework it is possible to discuss the effect of non-linear dispersion on the spectrum and 
correlations of (1). We add an extra term to the TL hamiltonian: 

H' = H + H^u (32) 

iJni = ^ / dxdx'h{x - x')yi;l{x)y^p{x'), (33) 
P ■' 

The subscript 'nl' stands for 'non-linear' dispersion. Function hq which is the Fourier transform of h{x) has these 
properties: 

/iq = Up < uf/A for \q\ < A, (34) 
hq < VF/\q\ for \q\ > A. (35) 

The second of these two conditions guarantees that for v'-p < v-p / A the modified kinetic energy pv^q + hqq'^ has the 
same sign as the original kinetic energy pvpq. That is, H^i does not induce an instability of the ground state by 
creating spurious Fermi point. 

Condition (34) implies that for small \q\ < A it is permissible to use 

Hni = v'^Y.I dxW^Pl{x)W%{x) (36) 

instead of (33). The quantity 2?;p is the Fermi velocity derivative with respect to the momentum. 

We need to express i?„i in terms oiij). To do so it is convenient to rewrite ifni with the help of density operators 
Pp rather than field operators V'p- This trick was already used in (19). It is easy to establish that 



4 {^^l{x)){V^p{x)): -\v^pp{x) = (37) 
p 

2ttv^ E 1™ {wpp{x) [:ipl{y){Vipp{y)):-: iVipl{y))ipp{y):] - 4Trbp{x - y)pp{y)} . 



P 

According to (19) the expression in the square brackets proportional to the product of two density operators. Thus: 



Hni(x) = 4 (VVi(a;))(VVp(a:)): -l^^Ppix) = (38) 
p 

V Km lim {pp{x) [pp{z)pp{y) - bp{z - y)] - 2bp{x - y)Pp{y)} ■ 

'-^ y^x z^y 



Now wc have to substitute ppq = uppq + vp-pq for g 7^ in the above expression. The resultant third order polynomial 
of p's must be rewritten in terms of tp with the help of (14) and (38). The final expression for H^i is: 

ffni = E / {v'f- {Vi>l)iVi>p): +ipg'p-p (: i>l{Vi>p): -: iVi>l)iJp: ) + ilpp} . (39) 
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The details of the derivation together with the exact formulas for coefficients , g' and fl can be found in Appendix 
C. Here we quote only the expressions which are valid if the interaction parameter is small. Let us define as 
Q!o = Q!g=o- When interaction is small (ao « go/^-KV-p <C 1) the coefBcients are: 

v'-p-v'-pK. v'pal, (40) 
ff'«27r4«o = ^, (41) 
il^-tK\'pal, (42) 

where 7 is a non-universal constant of order unity. 

The total hamiltonian H' is the sum of H, cq. (31), and i?ni- Due to iJni the quasiparticle dispersion becomes 
non-linear (the first term of (39)). Also, H^i introduces additional interactions between the quasiparticles. The 
operators corresponding to these interaction (the second term in (39)) are ©(wp^o) and irrelevant. Therefore, they 
can be neglected provided that 

g'h/v-p < 1 <^ v'pgoK < 4- (43) 

The quasiparticles remain free even if the bare dispersion is not linear as long as the bare interaction is sufficiently 
small. Note, that the free boson representation for TL model is much less tolerant of H^\. This term is equivalent to 
cubic interaction between bosons (see (38) and Ref.[l]). For the bosons to remain free a stricter condition WpA <C vp 
has to be satisfied. 

Equations (31) and (39) allow us to describe the thermodynamics of the electronic liquid with the generic dispersion. 
For example, the free energy equals to 

F = log (1 + e--^W/^) =^f log (1 + e--/T) ^ (44) 

where ep(fc) = pvpk + Wpfc^ is the dispersion of the quasiparticles. This result is non-pcrturbativc in WpA/fp. Cor- 
rections to the free energy due to the neglected quasiparticle interaction are ©((wp^o)^)- Thus, our result for F is 
accurate provided that (43) holds true. 

The specific heat at constant chemical potential can be derived from (44) in the limit of low temperature (T <C 



C{T) = ^T+ + O {{v'pT/vlf) + 0{gt) + 0{{g'n (45) 

We observe that non-zero dispersion curvature generates contribution to the specific heat. (Note: irrelevant 
operators neglected along the way also contribute to the specific heat. However, their contribution is 0{gQ) and 
0{{g'Y). For small go such corrections can be disregarded.) 

The same result for the specific heat could be obtained with the help of bosonization technique. There one must 
expand the free energy in orders of H^\. These calculations are done in Appendix D. They provide an important 
consistency check for the proposed approach. 



III. NON-RIGOROUS OVERVIEW OF THE METHOD 



We have finished the introduction of our method and ready to apply it to Green's function calculations. At this 
point we would like to make a break between two rather technical parts and explain why the method works on intuitive 
level. 

The central issue of this discussion could be loosely formulated as follows: why the strongly interacting bosons 
could be mapped on the weakly interacting fermions, what kind of 'cancellation' of interactions takes place? To 
answer this question we direct our attention to most essential aspects of the approach. All technical complications 
will be disregarded: we will ignore zero modes, normal ordering and assume that A = 00. This makes the following 
presentation more transparent. 

Let us start with the bosonized form of the TL hamiltonian, familiar to many researchers [1]: 

H' = Hkin + Hint+Hnl, (46) 

H^in + ifint = J dx^ {/C (Ve)' + /C-1 (V$)'} , (47) 
H^i = / dx^^ {(Ve + V^f + (VG - V$)'} , (48) 
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where boson fields © and $ are connected to the density operators in the usual way: 

Pp = -^(V$+pVe). (49) 

Tomonaga-Luttinger liquid parameter K. is defined below. For the purpose of this Section it is enough to remember 
that for small parameter K is close to unity: 

1-/C = 0(5o). (50) 
As we see from (48), the boson interaction is proportional to Wp. Now we perform the Bogoliubov transformation: 

Q = 1C^'^Q, (51) 
$ = /C-i/2$. (52) 

Our hamiltonian becomes: 

- \ 2 / _ \ 2 



We rewrite it as follows: 



^' = / '^^yjl^®) +(^*) I ^^^^ 



H' = Htrcc + ^H. (54) 

iff,ee = y'rfa;^|(ve)' + (v3.)'| + ^|(ve + v$)' + (ve-vl.)'|, (55) 

AH = J dx^^ I Ve + V$) ' + Ve - /C^/^ V$) ' (56) 

- (ve + v$)^- (ve- v$)^|. 

We split the above hamiltonian into two part for a reason. The first part -fffree becomes a free quasiparticle hamiltonian 
upon refermionization. Indeed, the first term in the curly brackets of (55) refermionizes to become the kinetic energy 
of the fermions with linear dispersion, the second term of (55) becomes q^-correction to the linear dispersion. 

Refermionized hamiltonian AH contains all the quasiparticle interaction. It also contains terms quadratic in tp 
which renormalize the value of v'p. The important point, however, is that AH is small if go is small: 

AH = 0((1 - IC)v'p) = O{gov^) (57) 

This can be established by observing that AH vanishes when go = ]C = 1. Note also, that AH is small if (43) 
is satisfied. Power counting shows that AH is irrelevant. Smallness and irrelevance of the quasiparticle interaction 
implies that the quasiparticlcs could be viewed as weakly interacting. 

We could look at our method from a different prospective. The hamiltonian H' = ifkin + Hmt + in addition to 
a quadratic (in if)) part ifkin has small marginal operator Hint. Plus, it has an irrelevant operator H^\ which we want 
to account for non-pcrturbativcly. Since perturbation theory in the marginal operator diverges, we must cither sum 
certain diagrams to all orders or, as it is done in the theory of superconductivity, perform a Bogoliubov transformation 
which kills the undesirable marginal operator. This is what our transformation U does: C/(i?kin + -ffint)t^^ is quadratic 
in fermionic fields. 

Operator _ffni is quadratic in ip. In general, however, operator UH^iU^ does not have to be quadratic. Thus, after 
the removal of the marginal operator by the Bogoliubov transformation, new interactions between the quasiparticlcs 
are generated. Fortunately, they are irrelevant and small. The smallness of the generated interactions becomes obvious 
if we observe that for small go transformation U is close to the identity transformation: U = 1 + O{go). Thus, the 
quasiparticle interaction induced by the action of U on H^\ could be treated perturbatively. This is the core of the 
approach we proposed in the previous Section. 
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IV. DENSITY-DENSITY PROPAGATOR 



In this section we will calculate density-density Green's function. The derivation of the correlation function for 
the total density operator p — pi^ + is quite simple if wc note that p is proportional to the quasiparticle density 
operator: pq = {ug + Vg)pg. Using this identity one obtains for small |g|: 

Vg{T) = -i (T. {p,(r)p_,(0)}) = /C (Pl,(t) + Vngir)) , (58) 

IC={ug+Vg) =W- ■ . (59) 

The chiral Green's function Vpg{T) = — {Tt {f>pq{T)pp-q{Q)}) /L will be calculated below to the zcroth order in the 
quasiparticle interaction constant g' . Unlike expansion in orders of Qq the perturbative expansion in orders of g' is 
a well-defined procedure:_ as we explained, the quasiparticle-quasiparticle interaction is irrelevant. With accuracy 
0{{g')'^) the propagator Dpg can be expressed as a convolution of two single-quasiparticle propagators: 

= "ZE Wi(fe+9)(0)}) {Tr , (60) 

k 

_ njipik)) - n(£p(fc + q)) 

ic. + £,(fc)-£yA: + g) • ^^'^ 

In eq. (61) n(e) is Fermi distribution function. Thus, at T = 0: 

- p ( iio - pvpq + v^q'^ \ 

= Iw:; ^,„-.,,,^.-.y n2 ' (62) 



A-Kv'-pq yibj — pvpq — v'^q'^ 
JC f iv' + iipq-v'^q^ ] 

This result has the accuracy of C((up5o)^), as explained above. The omitted corrections are due to the quasiparticle 
interaction g' . At the same time, the above expression for V is accurate to all orders in the dispersion curvature v'^. 
The retarded propagator Dg^ is obtained by the analytical continuation: 



For vanishing v'p we can write for the Green's function the following expansion: 

= + ^ - i.^,Lm' ) + ° V o («,.)^) . (65) 

The first term coincides with the well-known bosonization result [2] . The second term could be also found within the 
bosonization approach: one has to develop a perturbation theory expansion in powers of v'p. How this could be done 
is shown in Appendix E. These calculations serve as yet another consistency check for our method. 

One can make another interesting observation when examining (65). By looking at this expansion it is impossible 
to guess that the propagator (64) has a branch-cut, not a pole. To determine the correct complex structure the 
non-perturbative in Vp calculations are required. 

Such complex structure of the propagator indicates that there is no coherently propagating bosonic mode [2]. 
Instead, we have a continuum of quasiparticle - quasihole pair excitations. It is possible to visualize this continuum 
by calculating the spectral function Bgi^: 

Bg^ = -2lmDg^ = 7^{^ {^^ - {vFq - v'^q^f) - ^ {uj^ - {v^q + v^q^) } sgn uj. (66) 

(Note: unlike single- fermion spectral function Bg^^ does not have to be positive.) Every point of (g, w) plane where 
Bqu ^ corresponds to a quasiparticle - quasihole excitation with total energy and total momentum q. The set 
of these points is shown on fig. 2. 
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A previous attempt to account for non-zero Wp has been made in [1, 5]. However, the bosonic representation used 
in the latter references is not very convenient for such a task. To iUustrate the nature of the problem we consider a 
case of free electrons {go = 0) with strongly non-linear dispersion VpA ~ vp- 

Hfree = H^in + -f^nl- (67) 

The quasiparticle representation is trivial: U = 1 and ^ — ^. The spectrum of (67) is 

ek=VF\k\ + v'pk'^. (68) 

Density-density correlation function is given by (63). In the bosonization framework operator iJni corresponds to 
cubic interaction between bosons with the dimensionless coupling constant v'yA/vy of order unity (see (48)). To 
compute either spectrum or Green's functions one must resort to the perturbation theory whose accuracy, however, 
is not obvious due to lack of small parameter. 

If in addition to v'^ we have (?o 7^ then the Bogoliubov rotation of the TL bosons is required. When the Bogoliubov 
transformation acts on Hn\ it generates extra interaction terms of the form PppLp. The coupling constant for this 
kind of interaction is of the order of Wp^o- Thus, in the bosonic representation one has to deal with two kinds of 
interaction terms and two coupling constants one of which is of the order of unity. 

As we have seen in Section III, if we add weak electron-electron interaction Hint (eq.(3)) to H^ee (eq.(67)) the 
quasiparticle representation evolves continuously with (/o- U = 1 + O{go), H' = Hf^ee + C'(5o)- Non-zero go gen- 
erates interaction between the quasiparticles. Yet, this interaction remains small (if (43) is fulfilled) and irrelevant. 
Consequently, the Green's function is given by (63) with renormalized v-p and v'p, the spectrum of the interacting 
hamiltonian has the form (68). In short, the advantage of the quasiparticle representation steams from its ability to 
account for v'p non-perturbatively. 



V. SINGLE-ELECTRON GREEN'S FUNCTION 



Now we will show how to calculate the single-electron Green's function. At present we can find the Green's function 
for the case oi v'p =0 only. Although, this Green's function has been evaluated before by a number of approaches we 
would like to include this derivation to demonstrate different aspects of the method. 

In order to calculate the single-electron Green's function it is necessary to know how to express the electron field 
operator tp^ in terms of the quasiparticle operators and p. The following derivation answers this question. We 
begin by introducing unitary operator exp{(pp{x)) where 

(Pp{x) — ^^e^'^'^ ppq, f^pq — real coefficients. (69) 

Let us calculate the commutator of this operator with ppgi 

= (e*'''(="Vp9e-^''(^) - ppgj e'^''^^). (70) 
The first term in brackets can be calculated easily: 

e'^^'^^Vpge-^-^-) = ppg + [Mx),ppg] = ppg + pp-gpe-'o-. (71) 
Thus, the formula for the commutator is: 

[e'^-,Pp,]^pf3p^,e-'o-e^-. (72) 

If we choose Ppq = p the product tfjj,e'^'' commute with ppiq for any p' and g ^ 0. Therefore, it commutes with f2, 
eq.(9). This means that such product is invariant under the action of U. The action of U on the field operator is 
given by: 

= c/t^t[/ = ^j^le'^pUU-'^^U. (73) 

The latter formula is easy to invert: 

^l{x) = ^1,(:r)e'^-(")[/e-'^-(")i7t = i,l{x)Tl{x), (74) 



pVp(a:) „ 
^ 1 Fpq 
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= 0, (76) 

Wq=Uq- 1. (77) 

It gives us the desired equation for ip^ in terms of tp^ and p. 
Let us calculate the correlation function: 

(4(x,r)^p(0,0)) = (v;^(a;,r)^(x,T)^(0, 0)^^(0,0)). (78) 

The simplest way to evaluate this expression is to transform it in the following manner. The density operators with 
puq < 0{pnq > 0) must be shifted to the left (right) end, the quasiparticle field operators stay in the middle. The 
reason for such choice is that ppq |0) = for pUg > and (0| ppg = for pUg < where |0) is the ground state of (1). 
The details of this derivation are given in Appendix F. The result is 

r)%iO, 0)) = (^t .)4(o, 0)) cxp I - ^(1 - e^--^l^'^)^ } ^ J_ {^^) ^ ' ^''^ 

where wc used the notation 9 = 2w^|g^o = (^ + — 2)/2. The same formula is derived using bosonization [1, 2]. 

It is well-known fact that the above Green's function does not have pole. This means that the quasiparticle state 
has zero overlap with the physical electron state. 

VI. DISCUSSION 

In this paper we solve the TL model with the help of the unitary transformation. The transformation maps the 
original hamiltonian on the hamiltonian of weakly interacting quasiparticles. 

Our approach easily incorporates deviation of the electron dispersion from the linear form. As long as the bare 
interaction constant is sufiiciently small two new non-perturbative results (eq. (44) and eq. (63)) can be derived. The 
derived results accounts for the curvature parameter to all orders; they contain errors of order {g')^ ^ (vpgo)^ due 
to neglected interaction among the quasiparticles. 

The ability of our method to account for v'p to all orders allows us to resolve the complex structure of the density- 
density propagator. 

In principle, our diagonalization technique is not the only way to derive the quasiparticle representation. It is 
possible (fig.l) to bosonize (1) then perform the Bogoliubov transformation and then fermionize the diagonal bosonic 
hamiltonian [6, 7]. Alternatively, one can construct non-fermionic excitations using appropriate exponents of bosonic 
operators [8]. In our case use of bosonization looks like an unnecessary detour. 

To conclude, we propose a unitary transformation which maps TL model on a model of free fermions. Such approach 
reproduces or generalizes the TL correlation functions calculated using bosonization. 
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VIII. APPENDIX A 

In order to differentiate products of field operators one must perform normal ordering of the product under consid- 
eration. The subject of this Appendix is the definition of the normal ordering procedure we use in this paper. 

Most commonly a normal ordered product of two or more field operators is defined with the reference to some 
non-interacting ground state |0). For example, for product of two operators such definition reads: 




il;\x)^{x'y.= ^''{x)ip{x') - (0| ^''{x)ip{x') |0) 



(80) 
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The ground state is not specified in this equation. Thus, we actuaUy have infinite number of normal ordering 
procedures, each corresponding to a particular state |0). Usually, the problem at hand dictates the choice of this 
state. 

Normal ordered product (80) possesses two properties: (i) it is well-defined and analytical at a; = x'\ (ii) its ground 
state expectation value is zero. 

Since the ground state of the Tomonaga-Luttinger hamiltonian cannot be approximated by a ground state of non- 
interacting physical fermions TL ground state cannot be used in (80). Consequently, it is necessary to explain what 
state we use in our definition of the normal order. 

In the bosonization framework it is customary [1, 4] to normal order with respect to the the ground state of the 
hamiltonian; 



llQ = ^VYV}ic\c^. (81) 



fc 

We accept this definition in our work as well. Thus, for a product of two field operators we have: 

: ^l{x)^^{x'): = ^l{x)^^{x') - s^{x - x'), (82) 



V_ 

2'jTi{x — ipO) ' 



where Sp is equal to (0| ^pj,{x)^pp{x') |0). In this equation the pole of 'iljj,{x)'ipp{x') a,t x = x' is explicitly shown. The 
normal ordered product is well-behaved at x = x' and it satisfies the property (i), formulated above. 

It is not a miracle that property (i) remains intact despite the fact that the ground states of (1) and (81) are 
drastically different. Note that (i) is 'ultraviolet' in its nature - it refers to short distance (|a; — x'\ <C 1/A), high 
energy (|fc| 3> A) behavior of the operator product. Since interaction (3) is limited in g-space high-energy structure 
of the TL ground state is the same as that of (81). This guarantees that (i) is satisfied. 

Yet, the TL ground state expectation value of (82) for generic values of x and x' is not zero. That is, property (ii) 
is violated. Fortunately, we never need it in this paper. 

The normal ordered product of four field operators is defined by the equation: 

i,l{x)^p{x')^l{y)i,p{y') = : V^(x)Vp(x')^;(j/)^p(?/): +Spix - ?/): ^p(x')^;(2/): + (84) 

Sp{x' - y):ipl{x)ipp{y'): +Sp{x - y')sp{x' - y). (85) 

As above, the normal ordered product can be differentiated everywhere. 

For completeness, let us define the normal ordering of the density operators Pp(x): 

■ Pp{x)pp{y): = Pp{x)pp{y) -hp{x-y), (86) 
■■Pp{x)pp{y)pp{z): = pp{x)pp{y)pp{z)-hp{x-y)pp{z)-bp{x- z)pp{y)- (87) 

bp{y- z)pp{x), 

bp{x) = {sp{x)f. (88) 

We will use these definitions in Appendix C. 

The above rules can be generalized for T > 0: to normal order a product of density operators at finite temperature 
one has to replace function bp(x) by function: 



bpT{x) = {pp{x)pp{0)) = --^^^—j—y-j-^ (89) 



4{;| sinh"^(7rr(x/'DF — ipO)) 
It is obvious that bpx = 6p at T = 0. In Appendix D we will need the finite T definition. 

IX. APPENDIX B 

In this Appendix we provide the derivation of (30) starting from the formula (27). To express the hamiltonian (27) 
in terms of ^ rather than p we introduce a function A^: 

Aq = ul + v''^ + -^UgVg - 1 (90) 
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and its Fourier transform A(a;): 



(91) 



We defined gq such that gq ^ when \q\ oo. Thus, Vq ^ 0, Uq ^ 1 for \q\ ^ A. Therefore, the function A, 
vanishes for large \q\. This impUes that A(x) vanishes for |x| ^ 1/A. We may write: 



U HU'^ = ^ ^ PpqPp-q + ^ AqPpqPp. 



(92) 



pq 



pq 



ivp / dxy^p: jjlytjJp ■■ +~7^ / dxdx' K{x - x')pp{x)pp{x'). 

The first term in the above equation was obtained by inversion of (19). The second term must be normal ordered 
first: 



TTWf 



j dxdx'A{x — x')pp{x)pp{x') = — ^ X^ j dxdx' A{x — x') {bp{x — x')+: pp{x)pp{x'): } . 



(93) 



The additive constant proportional to / dxdx' Ah p will be disregarded. The normal ordered product of two p's can be 
expanded into Taylor series in powers of {x — x'): 



/ dxdx'A{x-x'):pp{x)pp{x')\ = ^ dx'A{x') / dx:pl{x): 
p p L J 

+ \ j dx'x''^A{x') j dx: Pp{x)V'^ pp{x): + . . . . 



(94) 



where ellipsis stand for higher order terms of the Taylor series. Since A(a;) vanishes for large |a;| the integrals of A are 
well-defined. We did not show the term proportional to -.ppVpp: since it is total derivative V:pp:, thus, it vanishes 

upon integration. 

We notice that the first term of expansion (94) is marginal. All other terms are irrelevant. Indeed, power counting 
shows that the scaling dimension of the operator : ppV^pp-. is 4 > 2. The scaling dimensions of the omitted terms are 
even higher. Thus, it is permissible to write: 



^0 = Aql^O- 



(95) 
(96) 



Consequently, transforming : Pp-. into the kinetic energy operator with the help of (19) and substituting the result into 
(92), we get the desired equation: 



UHU^ = /..^^p.p:^tv^^:+(irrelevant operators), 
vf = t;F(l + Ao) = t;F (^Uq + v; 



(97) 



.90 



(98) 



q=0 V V27ri'F, 

Note, that the omitted operators, in addition to being irrelevant, are also small (~ Aq ~ ffo/^p) interaction is small 
go < ^^F- 



X. APPENDIX C 



In this Appendix we show how i/ni can be expressed in terms of the quasiparticle field operators i^- We start 
with expression for the local density Wni: 



n^,{x) = ^;^5^:(VV^(a;))(VVp(a;)):-iv^Pp(a;) 



(99) 



Elim lim {pp{z) [pp{x)pp{y) - bp{x - y)] - 2bp{z - y)pp{y)} ■ 



3 "—^ z— »x 3/— ►x 

P 
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The last line ean be abbreviated if one use the definition of the normal ordering for the density operators pp{x), eq. 
(87). With this notation we can rewrite the formula for Ti^i: 

Wnl=^E^^r- (100) 

P 

To express the density operator Pp{x) in terms of the quasiparticle density operators Pp{x) we must remember that: 

ppg = UgPpq + VqP-pq for q 0, (101) 

Np = Np for zero modes. (102) 

Therefore, the density operator in co-ordinate space is equal to: 

Pp{x) = 6ppo + pI{x) + p%{x), (103) 

5ppO = ((1 - Uq)Np - VqN-^ , (104) 

where pp is the convolution of the density operator pp{x) with u{x): 

Pp{x) = {it* Pp){x) = I dx'u{x — x')pp{x'), (105) 

[uqe'"''. (106) 



l{x) = J^^ 



Likewise, p^p = ({* * p-p){x), where function v{x) is defined in the same manner as u{x). It is tempting to equate 
Hni and 



J2--i^ppo+pp+p-p)--- w 



p 



This, however, is not exactly accurate since the normal ordering and Bogoliubov transformation do not 'commute' with 
each other. Such phenomena becomes obvious if wc were to think in terms of Bosc creation and annihilation operators. 
Normal ordering places all creation operators to the left of all annihilation operators. Bogoliubov transformation mixes 
creation and annihilation operators, thus, it spoils normal ordering. 

To rewrite Hn\ in terms of p"'" correctly let us first examine the expression in square brackets in eq. (99): 

Pp{x)Pp{y) - bp{x - y) =: Pp{x)pp{y): = {5ppo + Ppix) + p%{x)) {Sppo + p^{y) + p"_p{y)) - bp{x - y) = (108) 
{: Pl{x)pp{y)- +{u *bp* - 2/) - Ki^ - y)] + {: P-p{x)p"_p{y): +(t) * 6-p * v){x - y)] + 
{pl{x)p%{y) + p%{x)pl{y)} + dppo {p^pix) + p%{x) + pl{y) + p\{y)} + Sp^^, 

where we applied the following transformations: 

Pp{x)Pp{y) = J dx'dy'u{x - x')u{v - y')pp{x')pp{y') (109) 

dx'dy'u{x - x')u{y - y') (: pp{x')pp{y'): +bp{x' - y')) 

=: pp{x)pp{y): +{u *bp* u). 

Similar transformations could be done for p"{x)p"{y). Since — l'^ = 1 the identity: 

u*f*u = f+[v*f*v) (110) 

holds true. Therefore: 

: Pp{x)Pp{yy. =: {p;{x) + p%{x)) {p;{y) + p%{y)) : + ({)* (&p + h.p) * {)) (x - y) + (111) 

(5ppo {~pI{x) + p%{x) + p;{y) + p%{y)} + Spl^, 



I' 
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The extra terms in this formula are non-singular functions of {x — y). In particular: 

{v*{b, + b.p)*v){x-y)^ J ^e^«(-^)|gb2. (112) 

Observe that the normal ordered product of p's differs from the normal ordered product of p's by a non-singular 
operator. 

Generalizing the above calculations for the product of three p's we obtain the following expression for HnV- 

H„i = — ^ ^: {p;f: +: {p%f: +3p^: {p%f: y6p\: {plf:+~c [dp^^ + ~p; + p\) + (z.m.), (113) 
p 

~c=^Jdq\q\vl (114) 

where 'z.m.' stand for zero modes terms which vanish in the thermodynamic limit. 

For small momenta the above formula can be simplified. What must be done, in its substance, amounts to inversion 
of equations (99) and (19). To illustrate this statement, consider the first and the second terms of (113). Limiting 
ourselves to the case A = oo we can neglect the g-dependence of Uq and Vq. Then, p^ = upp, Pp = vpp and it is 
possible to write: 

^ ^: {p; f: +: {p^f: = ^{u' + v') ^: p^: = v^u' + v') (: VV^tV^,: -^V^p,) . (115) 

p p ^ ^ 

When handling the third and the fourth terms of (113) we are to act in the similar fashion: the product : p^: should 

be rewritten in terms of ip with the help of (19). These two terms give the quasiparticle interaction. 

For a finite value of A our task becomes somewhat more complicated. In such a situation one can generalize the 
procedure of Appendix B. Let us briefly discuss the core of this generalization. As an example, consider the expression: 

: (pp3: = J dx'dx"dx"'v{x - x')v{x - x")v{x - x'"): Pp{x')pp{x")pp{x"'): . (116) 

It is easy to show with the help of Taylor expansion (see Appendix B) that the normal ordered product 
: pp{x')pp{x")pp{x"'): is equal to : Pp(a:;): +(more irrelevant operators), where 'more irrelevant operators' stands for 
operators whose scaling dimension is bigger than 3. Therefore, we obtain the following equation: 

• [ppf'- — ^' Pp- +(iiiore irrelevant operators), (117) 

where coefficient a is equal to f^|q=o- The expression : p^: has to be transformed further as in eq. (115). The remaining 
terms of (113) are transformed similarly. Therefore, one establishes: 

Wni = Y.h'^:{Vi,l)iV^p):-^V^Pp + ip~g'~p_p(^.,p^^^^ (118) 

P 

A(^L + Pr) + (more irrelevant operators), 

~g' = 27rv'p {ulvq+Uqvl)\^^^, (120) 

< (/ dq\q\vl^ . (121) 

The values of Uq and Vq are given by (23). Finally, we find for the hamiltonian: 

Hni = ^ / l^p: (Vi/'|)(VV'p): +ipg'p-p ^: V'p(Vi^p): -: {\7il)l)tpp: j + /ippj + (more irrelevant operators). (122) 
p •' 

We omitted the full derivative V^p from the hamiltonian since it contributes on the system boundary only. 

We don't want to prove here that the quasiparticle interaction produces no ultraviolet divergence of the perturbation 

theory. Instead, this issue is examined at the end of Appendix F. This is because the formalism developed there is 
particularly convenient for discussion of ultraviolet properties of the quasiparticle interaction. 
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Finally, let us briefly explain, why the quasiparticlc interaction generated by the action of U on H^i is small. This 
explanation is equivalent to the one given in Sec. III. Initially, our could be expressed as a cube of Pp{x), eq. (100). 
(Although, superficially, this expression appears to be sixth order in -0^, a subtler analysis shows that it is quadratic 
in ipp, eq. (99).) The Bogohubov transformation U converts (100) into (113). Unlike (100), eq. (113) is indeed a sixth 
order polynomial in ^pp. However, for small go operator U is close to the identity transformation: U = 1 + C'((7o)• 
Therefore, 

v'^:pl: = v:p:pl:+0{gov'^). (123) 

The first term on the right is proportional to Vipp^ipp, eq. (99). All generated terms, such as quasiparticle interactions 
and corrections to v'p, are O{gov'p). 

XI. APPENDIX D 

In this Appendix we show how the result for the specific heat, eq.(45), can be obtained within the bosoniza- 
tion framework. More specifically, 0{{vp)^) correction to the specific heat will be calculated with the help of the 
perturbation theory for the free boson hamiltonian: 

H=^^PpqPp-q- (124) 

pq 

We will see that it is the same as ©((wp)^) term in eq.(45). 

Our first step is to rewrite the non-linear dispersion hamiltonian H^i in a form particularly suitable for T > 
calculations: since in this Appendix we work at finite temperature it is natural to express Hni with the help of the 
finite temperature generalization of (86) and (87). As it was explained, such generalization is achieved by placing bpT, 
eq. (89), instead of bp, eq.(88). 

Using the new rules of the normal ordering we could prove that at T > the following is true: 

^fE^ Wp(^))Wp(^)):-^V%(a;) = l!^^:p3:+^p„ (125) 
p p ^ 

where : . . . : denotes here the finite temperature normal order. The proof of this equation is simple. Start from eq.(38): 

v'f {'7i>l{x)mp{^)): -l^'ppix) = (126) 

p 

^IJ^ V lim lim {pp{x) [pp{z)pp{y) - bp{z - y)] - 2bp{x - y)pp{y)} = 
6 '-^ y^x z^y 

p 

— V lim lim {pp{x) [pp{z)pp{y) - bpT{z - y)] - 2bpT{x - y)pp{y)} + 
6 '-^ y^x z^y 
p 

3pp{x) lim {bpT{x - y) - bp{x - y)} . 
y^x 

The last limit is 

, . r2 



lim {bpT{x -y)- bp{x - y)} = (127) 
Thus, bosonized Hni equals to 

H^, = ^ E / Pp- ) + ^(9ov'p). (128) 

This formula does not imply that Hni depends on the temperature. In this Appendix the symbol : . . . : denotes the 



finite T generalization of the normal ordering, consequently, the expression : p^: varies with temperature. The term 

compensates the temperature dependence of the normal ordered product : p^: so that H^i is temperature independent 
as it must be. We remind the reader that such a way of writing is a matter of convenience. 
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In the above expression we ignore the quasiparticlc interaction and the correction to the chemical potential since we 
are interested in ©((wp)^) corrections only while the neglected terms contribute to the thermodynamics at O{{vpgo)'^). 
The correction to the free energy due to i?ni is given by the usual perturbation theory formula: 

^^ = -\j dT{H^i{T)H^im = -ii'FfL j (iTdx(^i^(:pi(a;,r)::pl(0,0):) + ^(pL(a;,T)pL(0,0))) . (129) 

Using Wick's theorem we find: 

(: ~pI{x, t): : pt (0, 0): ) = 6(pl(x, t)pl(0, O))^, (130) 

(pL(a;, r)pL(0, 0)) = - . ^ — • (131) 

Av^ smh {■kT{x/v-p — it)) 

Therefore: 

f /tT^'T^ 1 TT^T^ 1 \ 

SF = (v-p)^L drdxi « « + «- 5 • (132) 

^ ^' J \ Q4 smh^{T^T{x/vF-iT)) SGv^smh'i'KTix/vF-iT))) ^ ' 



It is easy to calculate the integral 



/+00 
-00 



dx 2v-p 



sinh^(7rT(x/{;F - zr)) ttT 



for < r < /3. To find the integral of 1/ sinh^' wc use a trick. One can check directly that: 



sinh'' r 15 sinh^ r dr^ 
11 2 1 

20 sinh^ r 15 sinh^ r 



The integral of the full derivative is zero and we establish for < r < /3 

'■+°° -^^ e '■+°° dx IQVF 



r°° dx A / 

J-00 sinh^(7rT(a;/t;F - it)) ~ 15 J_, 



sinh^ {nT {x / Vf -ir)) 15 J-oo sinh^(7rT(a;/{;F - ir)) 157rr" 
Therefore, one finds the following corrections to the free energy, the entropy and the specific heat: 



^ djSS/L) _ UTT^{i'p)^T^ 
dT ~ 5i;| 

The last expression coincides with the second term of eq.(45) as it should be. 



(133) 



1 8 1 d'^Ur) 



^(^) = 7K-T4--T7-T2- (135) 



(136) 



77r^(f;^)^T^ 

30^' ^^^^^ 

,5/X = -^ = il^!(^, (138) 

' dT 15w| ' ^ ^ 



,C.r^ = i^^. (139) 



XII. APPENDIX E 

In this Appendix we calculate 0{{vp)'^) correction to the zero-temperature boson propagator. Such correction must 
coincide with the second term of (65). 

The correction in question to the chiral Matsubara propagator comes from ffni, eq.(128), with T = 0. The 
appropriate Feynman diagram is shown on fig. 3. As one can see from eq.(58) the corresponding expression is: 

5V^{x,t) = -| (^^) / dr'dT"dx'dx"{Tr {pp(0,0):p^(a;',r'): :pt(a=".r"):Pp(^,r)})„— + 0{{v'Fgo)^).{UO) 
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The time-ordered average with respect to hamiltonian (124) is equal to: 

J dr'dT"dx'dx"{Tr {ppiO,0): pl{x' ,T'y. : p^x" ,t"): Pp{x,T)}) co„- = (141) 

3x3x2x2 J dT'dT"dx'dx"V°{T',x') (v°{t" - t\x" - a;'))':D°(r - t",x- x"), 

where 'Dp(x, r) is the chiral free boson Green's function ■Dp(x,r) = — (l/Z/)(T-r{/Op(0,0)pp(x, r)}) for hamiltonian 
(124). 

The numerical factor on the right-hand side of (141) corresponds to 36 possible ways of contracting the density 
operators into a connected diagram presented on fig. 3. First, the external operator /5p(0, 0) could contract with either 
of two vertexes; the operator Pp{x, r) contracts with the remaining vertex. This gives a factor of two. Second, a given 
external operator could contract with either of three density operators in a vertex. This gives a factor of three for 
one external operator and another factor of three for the second external operator. Finally, after contracting with the 
external operators, every vertex has two uncontracted density operators. For them there are two ways of contraction. 
Thus, we have yet another factor of two. 

In the Fourier space we can write: 

-Dj,,^ « ICV%^ + 327r4/C(z}^)2 (p^J^^) ' 11^,^, (142) 

= ^ . , (143) 

27r iu> — pvpq 

Therefore, for correction to the full propagator we write: 



After analytical continuation we recover the 0{{vp) ) term of (65). 

XIII. APPENDIX F 

In this Appendix we provide a detailed derivation of the single-electron Green's function for TL model, equation 
(79). The formalism developed for the Green's function is also a convenient tool to prove that the quasiparticle 
perturbation theory has no ultraviolet divergences. 

We start with equation (78). First, the product J^^^p has to be normal-ordered: 

^\x,t)^ (0,0) = e"^p^'^+*™'°-'e~^p^^"'"'^^'°^e^'''-^^'''^'°'e~^p^^~*^'''°^ X (146) 
e-fE, ^r7(i-e""-'""')«',''>(P9)^pE, 7r7(i-<'""+'""')''?'»(-P9)^ 

Wp{x, x')=pJ2— fe^'"' - e''") w,ppg^i-pq), (147) 

Tin ^ ^ 

q ^ 

Yp{x,x') ^pJ2— (e*«^' - e^iA VgP-p.^ipq). (148) 

Tip ^ ' 

Here y = vpr. This equation was derived with the help of commutation relations (5). In order to calculate p(r) we 
used the expression (20). 

Second, we transform the whole expression V'j,(a^, t)^^{x, T)J^p{0, 0)tpp{0, 0). We shift exp(Wp) to the left past the 
quasiparticle field operator ip^. The exponent exp(— W|) is shifted to the very right. To perform these shifts the 
commutation rule (6) has to be used. 

ijl{x,T):Fl{x,r)^p{0,0)i,p{0,0) = {eW.(a.+iP.,o) (^t(a;, ^)^^(o, o)) e-wJC^+'fi/.o) } ^ (149) 
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After taking the expectation value of this expression and recalling that + 2wq = we obtain the desired expression 

for the correlation function. 

The above formula sets a convenient background for a rather general discussion of the ultraviolet properties of the 
quasiparticle-quasiparticle interactions. Imagine that we add the following term to the TL hamiltonian: 

^fextra = J2 [ ^xdx' Z{x - x')tPlix)'tJjp{x'). (150) 

p ■' 

Operator H-^x is a particular case of -f/oxtra with Z — —h". 

With the help of (149) it is easy to rewrite -ffextra in terms of the quasiparticle operators tp and p: 

H,^,,, = E / dxdx'Z{x - x')Cix - x') {eW.(-.-') (j;l{x)^^{x')^ q-'^IM | ^e^-^M^-^U-'-')^ , (151) 
P •' 

((x) =e^«^(^"''"'"'^'''. (152) 

Function ({x) is well-defined for |x| ^ 1/A and vanishes algebraically for large \x\. Operator (151) is normal-ordered, 
therefore, it is safe to expand the exponentials: 

H^^,.^= Y: H^mn'm'= E / ^^^^^ nlmfnMmM ^'^'^ ^^^^^ 



Since W = C(5o) ^^nd V = O(go) such expansion is controlled by the smallness of go. We want to argue that any 
interaction term Hnmn'm', n+m+n'+m' > does not lead to ultraviolet divergence of perturbation theory. Consider, 
for example, the second order correction to the ground state energy: 









/ Pl.-.PN 
1 \ fcl...fcjv 


HnOn'O 





k-i ...kjij 



Summation in the above formula goes over < n + r?,' + 1 quasiparticle momenta pi and over N quasihole momenta 
ki. We include m = m' = terms only because other terms annihilate the ground state. Thus, m, m' ^ terms do 
not contribute to the second order ground state correction. They do, however, contribute to higher order corrections 
and to corrections to different propagators. 

If the matrix element does not vanish at large momenta the expression (154) suffers from the ultraviolet divergence. 
We will show that the matrix element goes to zero if at least one \pi\ or |fcj| exceeds 2(n + n')A. 

Let us examine definitions of W and V. Since Vq and Wq vanishes for \q\ > A every individual operator W, V 
changes the total momentum by no more than A. Thus, momentum change produced by (W"V" ) is no more than 
(n + n')A. We will use this information in our next step. 

First step in evaluating the matrix element from (154) is creation of a quasiparticle-quasihole pair by applying tp^ip 
product to |0). Since the total momentum of all excitations is zero the momentum of this pair cancels the momentum 
of (W"V" ). Thus, the pair momentum magnitude is limited by (n + n')A. Two elementary excitations, quasiparticle 
and quasihole, of which the pair is composed have their momenta boimd by {n + n')A. 

As we explained above, by acting on the ground state the product ^p^^p creates an excited state with a quasiparticle 
and a quasihole. On this state the operator (W"V" ) acts creating yet another excited state. Every individual 
operator V or W when acting on a state with finite number of the elementary excitations can either (i) create 
another quasiparticle-quasihole pair or (ii) replace an elementary excitation with another one of the same chirality 
and charge but different momentum. 

In case (i) both newly created elementary excitations have their momenta bound by A. This is because single W 
or V operator cannot make quasiparticle-quasihole excitation with total momentum higher than A. 

In case (ii) an individual operator W or V can change momentum of an elementary excitation by no more that 
A. As above, this is a consequence of Vq and Wq vanishing at \q\ > A. Therefore, the product (W"V" ) cannot 
change momentum of an elementary excitation by more than (n -I- n')A. Consequently, momentum of an elementary 
excitation is boimd by 2{n + n')A. 

This means that in (154) summation is effectively performed over a sphere \ki\ < 2(n n')A, \pi\ < 2{n + n')A. 
Thus, this expression has no ultraviolet divergence. 

The above argument can be easily generalized to other perturbation theory formulas to prove that the perturbation 
theory for the quasiparticles has no ultraviolet divergences. 
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The absence of the ultraviolet divergences, together with irrelevance of the quasiparticle - quasiparticle interaction, 
implies that the quasiparticle perturbation theory is well-defined. 
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FIG. 1: Commutative diagram explains the relation between transformation U, and the bosonization. 
Bogoliubov transformation of bosons. 



'BT' stands for the 




FIG. 2: When = the spectral density of Dq^i is delta-function centered at a; = vpq line (dash line on the figure). This 
is the dispersion curve of the TL bosons. For Vp 7^ the spectral density is non-zero within the whole shaded area. This area 
represents the continuum of the quasiparticle-quasihole excitations. The TL bosons acquire finite life-time in such a situation. 




